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Deep Fabry-Perot Ha Observations of NGC 7793: 
a Very Extended Ha Disk and a Truly Declining Rotation Curve 

I. Dicaire^, C. Carignan^, P. Amram^, M. Marcelin^, J. Hlavacek-Larrondo^, 
M.-M. de Denus-Baillargeon^, O. Daigle^'^ and O. Hernandez^ 

ABSTRACT 

Deep Ha observations of the Sculptor Group galaxy NGC 7793 were ob- 
tained on the ESO 3.60m and the Marseille 36cm telescopes at La Silla, Chile. 
Ha emission is detected all the way to the edge of the H I disk, making of the 
H II disk of NGC 7793 one of the largest ever observed in a quiet non-AGN late- 
type system. Even in the very outer parts, the H II ionizing sources are probably 
mainly internal (massive stars in the disk) with an unlikely contribution from 
the extragalactic ionizing background. The Ha kinematics confirms what had 
already been seen with the H I observations: NGC 7793 has a truly declining 
rotation curve. However, the decline is not Keplerian and a dark halo is still 
needed to explain the rotation velocities in the outer parts. 

Subject headings: instrumentation: interferometers - techniques: radial veloci- 
ties - galaxies: kinematics and dynamics - galaxies: individual (NGC 7793) - 
galaxies: ISM 



1. Introduction 

NGC 7793 is one of the 5 bright members of the classical Sculptor Group along with 
NGC 55, NGC 247, NGC 253 and NGC 300. Sculptor is the closest group of galaxies outside 
the Local Group. It covers an a rea of ~ 20° in diamete r centered on the constellation Sculptor 
at a = 0'^30'" and 6 = -30° dde VaucouleursI Il959l ). However, it is still debated whether 
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Sculptor is really a group. iKarachentsev et al.l (120031 ) suggest instead that it may be part of 
a filament that extends along the line of sight from the Local Group out to ~ 5 Mpc. 

The optical parameters of NGC 7793 are summarized in Table [TJ It is a typical Sd 
galaxy with a very tiny bulge and a filamentary spiral structure (see figs 1 & 2). It is an 
intrinsically small syster n with an exponential sc ale length ^ 1 kpc. For this study, 
a distance of 3.38 Mpc (jPuche fc Carignanlll988l ) is adopte d, based on numerous dist ance 
indicators. However, sornewha t larger distances of e.g. 3.82 (IKarachentsev et al.l 120031 ) and 
3.91 Mpc ( jKarachentsevI l2005l ) have also been suggested. For our adopted distance, NGC 
7793 has an absolute B magnitude of -18.3 for a total blue luminosity of ~ 3.1 x 10^ . 

Previous ki nematical studies of NGC 7793 we re done in the optical, using Fabry-Perot in- 
terferometry by lDavoust &: de Vaucouleurj (Il980l ). Their Ha rotation curve (RC) extends to 
~ 4', which barely reaches the maximum velocity. Ten years later, deep VLA H I observations 
(13.5 ho urs in C/D configuration ) allowed one to derive the rotation curve twice as far, out 
to ~ 8' (ICarignan fc Puchdll990l ). A remarkable result of those H I observations was that, 
contrary to most spirals, the rotation curve is not fiat in the outer parts but appears to be 
declining {AVrot — 30 kms~^or 25% of Vmax between the maximum velocity and the last 
point of the RC), even after a careful modeling (tilted-ring model) of the warp (APA~ 20°) 
in the outer H I disk. While the RC is declining in the outer parts, it is less steep than a pure 
Keplerian decline and a dark halo is still needed to properly model the mass distribution. 

It would be interesting to be able to check the uniqueness of this RC with independent 
observations using another trace r than H I . As we have seen, previous Ha observations 
( jPavoust fc de Vaucouleurd Il980l ) only extend to ~ 4' or 80% of the optica l radius i?25- 
However, deep Ha imaging with a 68 A (FWHM) filter (IFerguson et al.lll996l ) showed that 
the diffuse ionized gas (DIG, also referred to as warm ionized medium, WIM) component in 
NGC 7793 easily extended out to i?25- Thus, one of the main motivations of the observations 
presented in this paper was the hope to extend the Ha kinematics, at least out to that point, 
in order to confirm (or not) the decline of the H I RC. 



The other motivation was to try to detect that DIG component out to very low levels, if 
possible at radii larg er than R^^. Large D IG components have been seen in nearby starburst 
and active galaxies (jVeilleux et al.l l2003l ) and can sometimes have sizes comparable to the 
H I component (see e.g. NGC 1068), but not much is known on the extent of the DIG 
component of quiet late-type systems. Comprehensive studies of the large scale structure of 



DIG components have be en done on edge-on systems (IRanc 



1992; Veilleux et al 



2006 



Oev et al.ll2007l ). 



1995h. on selected small regions of M31 dWalterbos fc Braunlll994h. and 



on la r ge samples (iKennicutt et al.lll995l : IZurita et al.ll2000l : iThilker et al 



et al.l 



199 



Dettmar fc Schulz 



200 



i 



Meurer et al. 
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Table 1: Optical parameters of NGC 7793. 

Parameter Value Reference 

Morphological type SA(s)d a 

R.A. (2000) 23^ 57™ 49^5 a 

Dec. (2000) -32° 35' 24" a 

Isophotal major diameter D25 10.1' b 

Holmberg radius Rho 6.1' b 

Exponential scale length a~^, kpc 1.1 b 

Axis ratio g = 6/a 0.60 b 

Inclination (go = 0.12), i 53.7° b 

Position angle, PA 279.3° b 

Corrected total B magnitude, B^^ {Ag = 0.02) 9.33 a,b 

Adopted distance (Mpc) 3.38 c 

(1' = 0.98 kpc) 

Absolute B mag nitude, M^'* -18.31 b 

Total blue luminosity {Mq = 5.43), Lq 3.1 x 10^ b 

(B-V) 0.54 a 



° lde Vaucouleurs et al.l (|l99ll ) 



Carignaij (|l985al ) 



iPuche fc CarignanI (|l988f ) 
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Massive stars clearly provide the largest source of Lyman continuum (Lyc) photons 
in spiral galaxies; an issue concerns whether the bulk of these photons are deposited over 
localized regions, such as the Stromgren spheres which define H II regions, or whether a 
significant fraction can escape from the regions of recent star formation w here they were 



created and thus ionize the interstellar medium over much larger scales (IFerguson et al. 



19961 ). The large scale radial distribution of the DIG across galactic disks can provide 
stringent constraints on the source of its ionization. For example, if the number of Lyc 
photons produced in star-forming regions were the only factor responsible for producing the 
DIG, then the ionized gas distribution should follow very closely, both on small and large 
scales, the distribution of discrete H II regions. This addresses the possible internal sources 
of ionization. 



However, if io nized gas is detected past the H I truncation edge (ICorbelli et al.l Il989 



van Gorkonj Il993l ). one possible explanation for that truncation could be that, at large 



ra dii, the thin H I disks get fully ionized by the metagalactic UV background, as discussed 



by iBland-Hawthorn et al.l (119971 ). This is another incentive to try to push the detection of 
the diffuse Ha as far out as possible. If Ha could be detected that far, it could even be used 
to put a limit on that UV background. 

Section II will present the new Fabry-Perot Ha observations and discuss data reduction 
while Section HI will derive the kinematical parameters and the rotation curves. This will 
be followed in Section IV by a study of the mass distribution and Section V will discuss 
the extent of the DIG. Finally, the principal results will be discussed in Section VI and a 
summary and the main conclusions will be given in Section VII. 



2. Observations &: Data Reduction 

The Fabry-Perot (FP) observations were obtained on the Marseille 36cm (from 2005 
October 27 to 2005 November 2) and on the ESO 3.6m ( 2005 Noveraber 3 ) telescopes at La 



Silla, Chile, using the photon counting camera FaNTOmM (iGach et al.ll2002l : [Hernandez et al. 



20031 ) . The same focal reducer, interference filter, Fabry-Perot etalon and camera were used 
on both telescopes. The interference filter, with a peak transmission Tmax ~ 80%, was 
centered at Ac = 6563A and had a FWHM of 30. 4A . The interference order of the Fabry- 
Perot interferometer was p = 765 at Ha for a free spectral range (FSR) of 8.6A or -390 
kms~^ . The mean Finesse of the etalon was ~ 18.5 (resolution ~ 14000). The finesse is a 

dimensionless value expressing the spectral resolution of the etalon (AA = fif§^)- The FSR 
was scanned in 60 channels by steps of 0.14A. The photon counting camera FaNTOmM is based 
on a Hamamatsu photocathode coupled with a Dalsa commercial CCD. The photocathode 
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has a quantum efficiency of ~23% at Ha and tfie CCD lias 1024x1024 12.5//m square pixels. 
Tlie CCD was operated in its low spatial resolution where pixels are binned 2x2 at 40 frames 
per second. 

A photon counting camera, such as FaNTOmM, is an essential tool for this kind of work. 
Its ability to rapidly scan the FP interferometer allows the photometric variations to be 
averaged out. For comparison, in CCD observations, each FP channel must be observed 
for at least 5 continuous minutes to avoid the read-out noise of the CCD from masking the 
weak galaxy's signal. This means that photometric conditions must not significantly change 
for ~5 hours with CCD observations. In photon counting, channels arc observed for 5 to 
15 seconds and a complete cycle is obtained every 5 to 15 minutes. Many cycles are made 
during an observation and since the data are analyzed on-line, SNR estimations can be 
made throughout the observations and the observer can decide when to stop the integration. 
The wavelength cahbration is done using a Ne lamp at 6598. 95A. Since the cahbration lamp 
is strong compared to galaxies' fluxes, calibrations are done in analogic mode (non photon 
counting). Typically, calibration channels are integrated 1 second each, such that a whole 
calibration only takes 1 minute. 

The galaxy was observed for a total of 100 minutes on the 3.60m (1.67 minutes/channel) 

and for 1200 minutes (20 minutes/channel) over 4 nights on the 36cm telescope. Raw observa- 
tional data consist of many data flies that contain photons' positions for every cycle/channel 
duo. With a cycling of 10 seconds integration time per channel, one file was created every 
10 seconds. The different steps of the data reduction are: 

• integration of raw data files into an interferogram data cube (3D data cube sliced for 
every Fabry- Perot channel); 

• phase correction of intcrferograms to create wavelength-sorted data cubes (3D data 
cube sliced for every wavelength interval); 

• banning spectral smoothing; 

• sky emission removal; 

• spatial binning/smoothing (e.g. adaptive binning); 

• radial velocity map extraction from emission line positions; 

• addition of astrometry information; 

• kinematical information extraction. 
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All the reduction was performed with IDL routines inspired of the ADHOCw reduction 
package (http://www.oamp.fr/adhoc/adhoc/adhocw.htm). Details on the data reduction 



can be found in iHernandez et al.l (120051 ) , IChemin et al.l (120061 ) and iDaigle et al.l (l2006al ) and 
the routines are available at http: / /www.astro.umontreal.ca/^odaigle/reduction. Especially, 
the way the profiles are found, the way the sky emission is subtracted, which is a crucial step 
when trying to detect very weak emission, the wa y the data is smoothe d using an adaptive 
smoothing technique are all described in details in iDaigle et al.l (l2006bl ). 




Fig. 1.— Fabry-Perot observations of NGC 7793 on the 3.60m. Top left : DSS Blue Band 
image. Top right : SPITZER IRAC 3.6/imimage. Middle left : H« monochromatic 
image. Middle right : Ha velocity field. Bottom : PV diagram. 
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Fig. 2.— Fabry-Perot observations of NGC 7793 on the 36cm. Top left : DSS Blue Band 
image. Top right : SPITZER IRAC 3.6/imimage. Middle left : Ha monochromatic 
image. Middle right : Ha velocity field. Bottom : PV diagram. 
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A first attempt at detecting the DIG in the outer disk of NGC 7793 was done in October 
1993, during 13 observing nights, for a cumulative total exposure time of 64 hours. The same 
instrumentation was used (36cm telescope, focal reducer, interference filter and interferom- 
eter), but the detector was different. This first data set was obtained wit h our old IPCS, 



a Th ompson camera based on a Silicium Nocticon TV tube technology (IBoulesteix et al. 



19841 ). However, the ca mera used to obta in the data set presented in this paper used a new 
GaAs tube technology (iGach et al.ll2002l ) and has a detective quantum efficiency ~25%, at 



least four times larger than the old Thompson IPCS. Thus, this second run was less time 
consuming (20 hours instead of 64) for reaching a higher SNR. Moreover, the shortening of 
the individual elementary exposure tim e together with the diminution of the total observ- 
ing time and better reduction software (jPaigle et al.ll2006bl ) allowed a better removal of the 
night sky lines which is critical at such a low signal level. 



3. Kinematical Parameters and Rotation Curves 

In order to derive the rotation curves, one must find the set of orientation parameters 
(rotation center {xq, yo), systemic velocity Vgys, inclination i, position angle of the major axis 
PA) that best represents the observed velocity fields, shown in fig. [1] and fig. [2], at all radii. 
The data in an opening angle of 40° about the minor axis are excluded from the analysis 
to minimize errors due to deprojection effects. A cosine weighting function (using the angle 
6 from the major axis) is used for the rest of the data, which gives ma ximum weighting on 



the m ajor axis. The task ROTCUR in the reduction package GIPSY (IVogelaar fc Terlouw 



20011 ) is used to find those parameters. 



Since they are correlated, the dynamical center and the systemic velocity are looked for 
first, by keeping i and PA fixed (using the photometric values given in Tabled]). We find that 
(xo, yo) corresponds to the optical center and that Vsys = 238 ±3 kms~^ for the 3.60m and 



230 ±4 km s^^ for the 36cm data. This can be compared to the ICarignan fc Puchd (Il990l ) 
H I global profile midpoint velocity of 230 ±2 kms~^ and intensity-weighted mean velocity 
of 235 ±4 kms~^ or to their systemic velocity derived in the same way from the H I velocity 
field of 227 ±3 kms^^ . The agreement is very good considering that a slight zero point shift 
can be expected since the calibration is done at the Ne wavelength and not at the observed 
wavelength. 

The next step is to obtain a least-squares solution for i, PA and Vrot in concentric annuli 
in the plane of the galaxy, by keeping (a;o,?/o) and Vsys fixed. This is shown in Fig. [31 For 
both data sets, an inclination i = 47° ±9° for the 3.6m and ±6° for the 36cm is adopted. 
For the PA we found 277° ±3° for the 3.60m and 286° ±4° for the 36 cm. This is consistent 
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with the H I data (ICarignan fc Puchdll990l ) which show that the PA is smaller in the inner 
parts than in the outer parts. In order to see if the solution found is a good representation of 
the whole galaxy, separate solutions are also obtained for the approaching and the receding 
sides. The 3 solutions (both sides, approaching, receding) and the adopted rotation curves 
are given in Fig. IHfor the 3.60m and in Fig. [5]for the 36cm data. The adopted errors are the 
biggest difference between the kinematical solution for both sides and the separate solutions 
for either the approaching or the receding sides or the intrinsic error from the tilted-ring 
model, if larger. The values of the RCs are tabulated in Table [2] and Table [31 respectively. 

For the 3.60m data of Fig. HI only the velocity points up to 132" will be used because 
further out there are not enough points on the receding side and the errors become too 
large. Similarly, for the 36cm data in Fig. 0, even if the Ha emission is detected out to ~8.5' 
(~ Rhi), only the data up to 7 ' will be used because there is no emission after that radius 
on the receding side. 



4. Mass Modeling 

As can be seen in Fig. |6l the Ha RC is very similar to the H I RC (errors on the 3.60m 
data were not plotted to facilitate the comparison; see Fig. S]) so that the mass models are 
expected to be very similar. However, since the H I velocities are slightly lower in the very 
inner parts and slightly higher in the intermediate parts, one could expect slight differences. 



In order to compare our results with those obtained using the H I RC (ICarignan fc Puche 



19901), the same approach for the study of the mass distribution will be used. The dark halo 
is represented by an isothermal sphere which can be described by two free parameters: the 
one-dimensional velocity dispersion a and the core radius Vc- The central density is then 
given by po = da"^ /AirGr'^. The other free parameter is the mass -to-light ratio (M /L)^ 



of the stellar disk. Details on the method used can be found in ICarignanI (jl985bl ) and 



Carignan fc FreemanI (119851 ) 



The parameters of the models and the main results are given in Table H] and illustrated 
in Fig. [3 As expected, the slightly higher Ha velocity ~1 kpc gives a higher (A^/L^)^ for 
the stellar disk, thus reducing the dark component and the slightly lower Ha velocities ~5 
kpc also gives a lower density dark component. Those two factors, combined with the fact 
that the H I RC extends slightly further out, explain the fact that the Mdark/Mium goes down 
from 0.97 to 0.75 at the last point. Nevertheless, within the uncertainties, one can say that 
the Ha data confirms the results obtained with the H I data with Mdark/Mium — 0.75 — 1.0 
at the last velocity point for a total mass for NGC 7793 of ~ 1.3 — 1.5 x IO^^Mq. Again, 
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Fig. 3. — Kinematical parameters of NGC 7793 derived from the 3.60m (left) and the 36cm 
(right) Fabry-Perot observations. 
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Fig. 4. — Left : Rotation curve derived from the 3.60m data for the whole galaxy (black) and 
separately for the receding (red) and the approaching (blue) sides with Vgys = 238 kms~^ , 
PA= 277° and i = 47°. Right : 3.60m adopted Ha rotation curve for NGC 7793. The data 
are also tabulated in Table [2l 
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Fig. 5. — Left : Rotation curve derived from the 36cm data for the whole galaxy (black) and 
separately for the receding (red) and the approaching (blue) sides with Vsys = 230 kms^^ , 
PA= 286° and i = 47°. The H I data are plotted as crosses (green). Right : 36cm adopted 
Ha rotation curve for NGC 7793. The data are also tabulated in Table |3l 
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Table 2: 


Ha RC from the 3.60m data. 






Kadius 


'rot 


'rot 


J- L;C*J\J.J. HO 


T r 


ATT" 


(arcsec) 


1 JY±1± ij I 


r km "1 


I CLL v^ov;v^ 1 


/I — 1 \ 

( km s ) 


/I — 1 \ 
( km s ) 


6.0 


Q 1 


8 5 


70 


67.3 


10.1 


8.0 


1 7 fi 


7 4 


79 


70.3 


9.7 


10.0 


IQ 3 


7 4 


74 

I t:. VJ 


69.6 


9.8 


12.0 


24 8 


6 4 


76 

1 VJ. VJ 


70.5 


9.3 


14.0 


30 


7 2 


78 


69.4 


8.8 


16.0 




7 


80 
ou. u 


71.3 


11.2 


18.0 


38 3 


7 6 


82 

kJ^. VJ 


73.2 


11.2 


20.0 


42 7 


8 8 


84 

VJ 


75.4 


8.8 


22.0 


47 


5 3 


86 

(JVJ. VJ 


76.1 


7.2 


24.0 


4Q 8 


5 6 


88 

VJ 


75.3 


11.4 


26.0 


52 2 


fi 2 


90 

(JVJ. VJ 


79.0 


11.9 


28.0 


53 4 


6 3 


92 

<J^j . VJ 


80.6 


11.1 


30.0 


54 fi 


5 Q 


94 

iJrt. VJ 


80.3 


13.6 


32.0 


55 5 


6 Q 


96 


80.9 


14.1 


34.0 


55 Q 


7 5 

1 . tJ 


98 

<JO . v^ 


84.6 


13.7 


36.0 


58 3 


7 4 


100 


85.5 


13.5 


38.0 


5Q 3 

\JtJ.KJ 


6 3 


1 02 

i V.7Zj .VJ 


88.1 


10.2 


40.0 


5Q 5 


6 7 


104 

i VJ'-± .VJ 


87.3 


10.5 


42.0 


61 1 


6 1 


106 

i VJVJ .VJ 


85.3 


9.3 


44.0 


61 8 


6 6 

VJ . VJ 


1 08 

± VJO . VJ 


86.7 


8.1 


46.0 


61 7 


8 


110 


84.4 


10.0 


48.0 


62 7 


8 5 

Kj . tj 


112.0 


83.7 


-1/1/1 

14.4 


50.0 


63.2 


8.1 


114.0 


89.9 


12.7 


oz.U 


62.4 


9.2 


116.0 


on Q 

yu.o 


iU.o 


54.0 


63.9 


7.4 


118.0 


89.9 


10.9 


56.0 


61.8 


9.3 


120.0 


92.0 


13.6 


58.0 


61.1 


11.4 


122.0 


91.8 


22.1 


60.0 


62.4 


10.2 


124.0 


94.2 


23.9 


62.0 


66.4 


9.2 


126.0 


93.0 


22.3 


64.0 


66.7 


10.7 


128.0 


95.4 


22.3 


66.0 


66.8 


10.2 


130.0 


94.7 


21.7 


68.0 


65.7 


10.6 


132.0 


87.9 


13.1 



-14- 



Table 3: Ha RC from the 36cm data. 
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Fig. 6. — H I (green), 3.60m Ha (blue) and 36cm Ha (red) rotations curves of NGC 7793. 
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one can see that even if the RC is dechning over the second half of the radius range, a dark 
halo is still needed to reproduce the observed kinematics. 



Ha Extent 



The first mention that NGC 7793 had an extended ionized component was by iMonnet 



(I1971I ) from quit e shallow photographic data. A deeper and more detailed study was done 
more recently by iFerguson et al.l (Il996l ). They showed that the observed Ha luminosity of 
the DIG component of NGC 7793 was large, accounting for 40% of the to tal Ha emission. 



simil ar to what is see n in M 31 (IWalterbos fc BraunI 1 19941 ) or the SMC (IKennicutt et al. 



I995I ). iFerguson et al.l (119961 ) argue that the integrated minimum Lyman continuum (Lyc) 
power required to sustain the DIG is enormous and can only easily be met by the ionizing 
output from massive star formation; the mechanical energy from supernovae and stellar winds 
falls short of that required by at least a factor of 2. Their results support the hypothesis 
that the DIG is photoionized by Lyc photons leaking out of discrete H II regions. 

A comparison between our deep Fabry-Perot Ha map and the deep Ha continuum sub- 
tracted image from Ferguson et al. shows a very good general agreement. Nevertheless, we 
reach fainter diffuse Ha emission in the outskirts of the galaxy locate d between H II r egions . 
In addition, spurs are visible on our image that were not detected by IFerguson et al.l (119961 ). 
Fig. E shows the Ha radia l profile of NGC 7793, derived using the GIPSY task ELLINT 
( IVogelaar fc Terlouwll200ll ). Since our data were not photometrica lly calibrated, our pro file 
was adjusted in order to be in rough agreement with the profile of IFerguson et al.l (1l996l ) in 
the outer parts. Th i s shou ld be looked at as just an order of magnitude representation since 
the IFerguson et al.l ( 19961 ) data were obtained in imagery with a 68 A filter which included 
not only Ha but also th e [N II ] lines. While th e y cor rected for the [N II ] contamination 
using [N II ]/Ha = 0.22, iBland-Hawthorn et al.l (11997 1) found a ratio close to unity in the 
outer parts of the other Sculptor Group galaxy NGC 253. So, even if the intensity scale 
could be off by a factor of a few, this is sufficient to show that our long Fabry-Perot ob- 
servations (20 hours total or 20 minutes per channel) on a 36 cm telescope nearly reach 



the limit of the H I disk {Dhi = I.7D25 ; ICarignan fc Puchelll990l ). To be a ble to reach so 



large radii is interesting since it was estimated that if the truncation edge (ICorbelli et al. 



19891 : Ivan GorkomI Il993l ) of the H I disks i s due to ionization of the hydrogen by the cosmic 
background radiation, the expected flux (IBland-Hawthorn et al.l Il994j . Il997l ) should be an 
Ha emission measure £m(Ha ) ~ 0.2 — 2.0 pc cm~^, which is the level reached by our deep 
observations. 



Contrary to what was found by 



Bland-Hawthorn et al, 



(I1997I ) (hereafter referred to as 
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Table 4: Parameters and results for the mass models. 



H I rotation curve" 



36cm Ha rotation curve 



Luminous component (yW/L^)^ = 2.2 A^q/L^ 







{M/Lb), = 2.6 ±0.3 A^o/Lq 



Dark halo component 



rc 
a - 

Po 



-- 2.7 kpc 
40.8 kms"^ 
= 0.038 MqPC-^ 



Tc = 2.9 ±0.8 kpc 

a = 37.0 ±6.0kms~i 

Po = 0.027 ±0.009 MqPC-^ 



At Rho (6.0 kpc) 



Phalo = 0.0025 MqPC 
Mdark/Mi^^ = 0.80 



M^aark+lum) = 1.3 X 10^° Mq M^dark+lum) = 1.2 X lO^^^ 



Phalo = 0.0020 A^opc- 

Mdark/Mi^^ = 0.67 



10 







At the last point: 



f outer- — 7.35 kpc 

phalo = 0.0015 A^opc-^ 

Mdark/Mi^^ = 0.97 



M, 



(dark+lum) 



1.5 X 10^" M 



router = 6.84 kpc 

phalo = 0.0015 MqPC-^ 

Mdark/Miurn = 0.75 



M, 



{dark+lum) 



1.3 X 10^" M 



" ICarignan fc Puchj |l990l ) 

^this work 
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2 4 6 8 

Radius (kpc) 



Fig. 7. — Mass model (isothermal sphere halo) for the 36cm Ha rotation curve (filled circles) 
of NGC 7793. The H I curve (open circles) is also shown. 
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BFQ1997) for NGC 253, the Ha disk of NGC 7793 does not extend past the H I emission 
even if it gets close to it. Indeed, BFQ1997 have detected ionized gas out to 1.4 R2^, while 
the H I extends only out to 1.2 i?25- We should nevertheless notice that these authors did 
not give an Ha map of NGC 253 but only presented some discrete Ha emission located at 
a single position. Using Fabry-Perot techniques as well, they computed a radial velocity 
for this emission region and since they used a fixed etalon, instead of a scanning one as we 
did, they do not have full spatial coverage of the galaxy. Moreover, they have only observed 
one side of the galaxy. Further observations of NGC 253 are needed to verify if the outer 
H II region detected by BFQ1997 is isolated or if it is a fragment of a more extended region. 
In addition, it is necessary to search for a counterpart of this region on the other side of the 
galaxy to strengthen their conclusions. 



Discussion 



The outer regions of galaxies are expected to deliver important clues on the way galaxies 
are assembled. It is into these regions that material was accreted during violent interactions 
in the past and where new material continues to fall today. Cosmological simulations of 
disk galaxy formation predict a wealth of baryonic structures and substructures up to ultra- 
faint surface brightness levels, possibly due for instance to tida l debris from accretion an d 
disruption of an expected rich population of luminous satellites ( iBuUock fc Johnstonll2005l ). 
The faintest features visible in M31 reach effective V-band surface brightnesses greater than 

~ 30 mag arcsec"^. The currently most powerful technique to probe the very low 
surface brightness regions is resolved star counts, but it is only applicable to our nearest 
large neighbors, such as M31 and M33. Both galaxies appear to show evidence for metal- 
poor, pressure-supp orted stellar halos but only M31 shows evidence for recent accretion 
( jlbata et al. 1120051 ). The comb i natio n of diffuse light surface photometry with resolved 
star counts allowed llrwin et al.l (120051 ) to reach an effective V-band surface brightness of 
Hv ~ 32 mag arcsec"^ for M31. Globular cluster and planetary nebula populations also offer 
additional clues to the formation and evolution at large radii. 

The presence of the neutral gas component, detected through its H I emission, at radii 
several times larger than the optical disk (measured from its surface brightness at /i^ = 25 
mag arcsec"^) has now been well known for more than 30 years. Many of the extended 
H I rotation curves are rising or flat at their outermost points, while others are slightly de- 
creasing. Nevertheless, none of them sho w a Keplerian decline, m eaning that the edge of the 
mass distribution has not yet been found (jPuche fc Carignanlll99ll ). The only suggestion that 
such an edge could have been reached is for the dwarf galaxy DDO 154 (jCarignan &: Purton 
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120 




Radius (R/R ) 



Fi g. 8. — Radial Ha pro file of NGC 7793. Superposed on it is the radial profile (filled circles) 
of iFerguson et al.l (Il996l ). 



- 21 - 



19981 ). Deep H I observations have shown that H I disks are abruptly truncated at column 
densities near 10^^ — 10^^ atoms cm~^, where the rotation curves still do not have a Keplerian 
behavior (e.g. Van Gorkom 1993). If the metagalactic UV backgroun d is sufficiently strong, 
beyon d a certain galactic radius the cold gas may become ionized ( iBochkarev fc Siuniaev 
19m 



The total Ha emission from star forming galaxies i s divided roughly in half, between 
discrete H II regions and the DIG. Pioneering work by iReynoldd (119841 ) established that 
only the massive OB stars are able to ionize the DIG. On the other hand, all the observers 
agree on the fact that the DIG fraction seems to be independent of galaxy parameters 
such as star formatio n rat e and Hubble type. Th is looks a bit controversial. More recently, 
Reynolds et al.l (119991 ) and I Collins fc Rand J200lh suspected that OB stars cannot exclusively 
explain the ionization state of the DIG. iDopita et al.l (120061 ) suggested that part of the DIG 
may consist of extremely evolved filamentary H II regions that may be difficult to detect. 

On the other hand, if ionizing radiation escapes from H II regions to ionize the DIG, it 
does escape from the galaxies and affects the ionization state of the intergalactic medium 
(IGM) and the galactic medium of external galaxies which may be partially at the origin of 
this UV radiation. Indeed, a key question is the origin of the ultra-violet photons responsible 
for ionizing the hydrogen in the intergalactic medium (the Lya forest). The most likely 
sourc es of ionizing photons are young hot stars in galaxies and QSOs (IBarkana fc Loeb 
20011 ) but their exact origin remains one of the unsolved questions of modern cosmology. 
Observations of high-r edshift quasars imply that the reionization of cosmic hydrogen was 
completed by z~6 (e.g. Fan et aul2006 . and reference therein). Lack of observational evidence 
and theoretical uncertainties makes it difficult to estimate the contribution made by both 
sources and leads to inconsistent results. To complicate the scheme, th e relative contrib ution 
of QSOs and star forming galaxies depends on the lookback time ( iFan et al.l |2006| . and 
reference therein). 

Recently, based on the first data release of the SINGG survey of Hl-selected galaxies 



( iMeurer et al.ll2006l ). lOey et al.l (120071 ) have undertaken a quantification of the ionization 
in neutral galaxies and the diffuse and warm ionized medium (NGC 7793 is unfortunately 
not included in the first SINGG data release of 109 galaxies). They concluded that the 



Thilker et al. 


2002: 


Zurita et al. 


2000 



Wang et al.l 



1999 



Hoopes et al.l Il996l ) and also consistent with the va l ue ~ 0.50 found by 



Ferguson et al.l (119961 ) for NGC 7793. The difference between lOey et al.l (120071 ) and previous 
studies may be due to selection effects favoring large, optically bright and nearby galaxies 
with high star formation rates. NGC 7793 is closer than any (except one) of the first 
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SINGG data release. We have computed, for NGC 7793, its H a effective ( half l ight) radius 
~ 1.97±0.05 kpc (~ 120 arcsec). Refering to equation (1) fromlOev et al.l (120071) . and using 



the Ha total luminosity of 4.1 10^° erg.s ^ computed by lFerguson et al. 
the logarithm of effective surface brightness hogiT^Ha] 



(Il996l ) , we found that 
38.65 Log(erg s~^ kpc~^. Following 
Oey et al.l (l2007l ). this brings NGC 7793 to be classified as a normal galaxy with respect to 



its star formation rate. 



6.1. Possible Origins of the Ionization: Internal and/or External 

Several mechanisms for the ionization of the neutral hydrogen and for heating the elec- 
tron population belonging to the DIG may be found in the literature. BFQ1997 review them 
in the specific context of NGG 253. The main source of ionization and heating may be 
(see references in BFQ1997): the metagalactic ionizing background; compact halo sources; 
compact disk sources (white dwarf populations or horizontal-branch stars, young enough to 
produce significant UV emission); ram pressure heating (NGC 253 and NGC 7793 lie both in 
the Sculptor group having little intracluster medium); turbulence-driven MHD-wave heating; 
mixing layers driven by bulk flows; galactic fountain material driven by the inner starburst 
through shocks; dilute photoionization; gas phase depletion (induced by the formation of 
grains); low energy cosmic-ray electrons and young stellar disk. 

BFQ1997 ruled out most of these mechanisms, except the ionization by the hot young 
stars of the disk. Even if those young stars are in the disk and most of them further in, 
the warp of the H I disk (which is the case for both NGC 253 and NGC 7793) makes it 
possible to expose it to their radiation. In particular, they ruled out the ionization by the 
cosmic UV background which is suspected to be ~5 times below the ionizing flux necessary 
to explain the emission seen in NGC 253. As the Ha disk does not reach the outskirts of 
the H I disk of NGC 7793, it is very unlikely that the source of ionization of the DIG is the 
metagalactic UV background. Indeed, if the metagalactic UV background was at the origin 
of the ionization it should have principally affected the outer H I regions, which is not the 
case. An alternative could be that the Ha detection level is not yet low enough to be able 
to detect the ionization by the cosmic UV background radiation, thus, this should mean 
that the expected flux should have an Ha emission measure £^m(Ha ) lower than what was 
measured for NGC 7793, namely ~ 0.2 — 2.0 pc cm~^. 
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6.2. A Truly Declining Rotation Curve 



The existence of falling RCs was first claimed in the co ntext of t he ef fect of cluster 
environments on the member galaxies. Following the work of iDresslerl ( 119801 ). the increase 
of intergalactic gas and of gravitational interactions expected toward the centers of clusters 
should not only affect t he morphologica l type s of the systems but also their kinematics 
and mass distribution. IWhitmore et al.l (119881 ). from long-slit emission-line observations, 
presented evidence for a correlation in the sense that the inner-cluster spirals tend to have 
falling RCs, while those of the outer clu ster and the v a st ma jority of field spirals seem to 
have fiat or even rising RCs. However, lAmram et al.l (119931 ). using full 3-D Fabry-Perot 
data did not confirm t hese correlations and moreover found that many galaxies for which 
Whitmore et al.l (119881 ) claimed a falling RC, had in fact a fiat or even rising RC after a 
careful analysis of the Fabry-Perot 2D velocity fields (lAmram et al.lll996l ). 



Declining RCs haye been seen before in other systems ( Casertano fc van GorkomI Il991 



Honma fc Sofuelll997l : lRyder et al.lll998l : ISofue &: Rubinll200ll : lNoordermeerll2007l ). However, 
those galaxies are most of the time earlier type systems (SO - Sbc) with highly-concentrated 
stellar light distribution and the maximum rotation velocity is reached at small radii (few 
hundreds kpc and < a^^). This is certainly not the case for NGC 7793 which has a very 
small bulge and a relatively diffuse (exponential) light distribution and where Vmax is reached 
at ~ 4a~^. Keplerian declin es have also been o bserv ed in a few ed ge-on systems ( e .g. N GC 
891: ISancisi fc All^ Jl979h : lOosterloo et"aD J2007h . NGC 3079: Iveilleux et all Jl995h & 
NGC 4244: lOUingl (jl996l )). but the interpretation in the case of those galaxies is much more 
difficult. Finally, a truncation of the halo at ~ 10 kpc was advocated to exp lain the decline 
of the last two points of the RCs of NGC 5204 Jsicotte fc CarignanI Il997t ) and NGC 253 
(IBland-Hawthorn et al.l 119971 ). but here the decrease of Vmax and the last velocity point is 



only ~ 10% of Vmax and ~ 10% of Vmax, while, for NGC 7793, the decrease is ~ 25% of Vmax 
and over about half the radius range. 

So, for late-type nearly bulgeless galaxies, the decreasing RC of NGC 7793 is rather 
unique. As can be seen in Fig. El the decline of the RC seen with the H I data is confirmed 
by the Ha 36cm data and the kinematical analysis of the 3 sets of data (H I & Ha ) give very 
similar results. Since a bulge component can not be invoked to explain the falling gradient of 
the RC, NGC 7793 really seems to have a genuine declining RC. The new data also confirm 
that the decline is slower than Keplerian and that a dark component is still necessary to 
explain the observed kinematics. 
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7. Summary and Conclusions 

Deep Ha observations of NGC 7793 liave been presented. Tlie main results are: 

• A total of 20 hours of observations was obtained on the Marseille 36cm telescope which 
allowed to reach sensitivities of the order of £m(Ha ) ~ 0.2 — 2.0 pc cm~^; 

• this allowed the detection of diffuse Ha emission out to the edge of the H I disk; 

• a warp of the disk plane, already seen in the H I data is also seen in the H II disk; 

• a RC is derived out to 2.2' from the 3.6m data and out to 7' from the 36cm data; 

• the two Ha RCs superposed exactly with the H I RC; 

• the Ha results confirm the H I findings that NGC 7793 has a truly declining rotation 
curve; 

• a model of the mass distribution gives a total mass for NGC 7793 of M(^dark+ium) = 
1.3 - 1.5 X 10^° A^ofor Mdark/Mi^rn = 0.75-1.0 at the last velocity point. 



This was our first experiment of doing such long integrations using a 36cm telescope 
coupled with our FaNTOmM scanning Fabry- Perot system. This yielded very interesting results 
for NGC 7793 by detecting diffuse Ha emission all the way to the edge of the H I disk and 
confirming its peculiar kinematics. In the near future, we intend to carry out the same kind 
of observations for all the Sculptor Group galaxies, hoping to reach even fainter levels . This 



should be possible with an even more sensitive camera soon available (jPaigle et al.l 12004 



2006d ). based on an EMCCD, which should allow us to increase the detective quantum 



efficiency by a factor ~ 3-4. 
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